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Abstract. Let T be a general complex tensor of format (ni, ...,n<j). When the fraction 
IT, — 1)] is an integer, and a natural inequality (called balancedness) is satisfied, 

it is expected that T has finitely many minimal decomposition as a sum of decomposable 
tensors. We show how homotopy techniques allow us to find all the decompositions of T, 
starting from a given one. Computationally, this gives a guess regarding the total number of 
such decompositions. This guess matches exactly with all cases previously known, and pre¬ 
dicts several unknown cases. Some surprising experiments yielded two new cases of generic 
identifiability: formats (3,4, 5) and (2, 2, 2, 3) which have a unique decomposition as the sum 
of 6 and 4 decomposable tensors, respectively. We conjecture that these two cases together 
with the classically known matrix pencils are the only cases where generic identifiability 
holds, i.e., the only identifiable cases. Building on the computational experiments, we use 
algebraic geometry to prove these two new cases are indeed generically identifiable. 


1. Introduction 


Tensor decomposition is an active field of research, with many applications (see, e.g., [49] 
for a broad overview). A tensor T of format (ni,..., n^) is an element of the tensor space 
C ni <E> ■ ■ ■ <8> C nd . The rank of T is the minimum r such that 


(1) T = J2 V 

i —1 

where v'- G C nj . This reduces to the usual matrix rank when d = 2. 

The space C ni <g)- • -<S>C nd contains a dense subset where the rank is constant. This constant 
is called the generic rank for tensors of format (ni,..., n^). By a simple dimensional count, 
the generic rank for tensors of format (ni,.rz.^) is at least 


( 2 ) 


R(ri \,... ,n d ) := 


Utim 


n 


i= 1 


rii 


1 + J2i=i( n i — i) i-^ + 52 = i 


71 ; 


The value \R(ni ,..., n^)] is called the expected generic rank for (77!,..., rid)- 

A necessary condition for a general tensor T of format (77!,... ,rid) to have only finitely 
many decompositions ([I]) is that the number R(ni,... ,nd) is actually an integer. Such 
formats are called perfect [13168]. Moreover, if a general tensor is known to have finitely many 
decompositions ([I]), then the generic rank is equal to the expected generic rank _R(t7i, ..., rid)- 
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A tensor format is said to be generically identifiable if a generic tensor of that format 
has a unique decomposition (up to global rescaling and permutation of the factors). A 
tensor format is said to satisfy perfect identifiability if the format is perfect and generically 
identifiable. The main goal of this paper is to study the number of decompositions of perfect 
formats (ni,... ,rid) when the generic rank is indeed equal to the generic expected rank. 

The main tool for inquiry is numerical algebraic geometry, a collection of algorithms 
to numerically compute and manipulate solutions sets of polynomial systems. Numerical 
algebraic geometry, named in |65] , grew out of numerical continuation methods for finding 
all isolated solutions of polynomial systems. For a development and history of the area, 
see the monographs [TUI 155] and the survey 133- The monograph uni develops the subject 
using the software package Bert ini [Sj, which is used to perform the computations in this 
article. For understanding the relation between numerical approaches and the more classical 
symbolic approaches to computational algebraic geometry, see |7]. 

Numerical algebraic geometry has proven useful in many other applications. A small subset 
of such applications include computing the initial cases for equations of an infinite family of 
Segre-Grassmann hypersurfaces in [24] : numerically decomposing a variety in HU which was 
a crucial computation leading to a set-theoretic solution of the so-called salmon problem [2] 
improving upon a previous result of Friedland m, and inspiring a later result of Friedland 
and Gross PH]; solving Alt’s problem [5] which counts the number of distinct four-bar linkages 
whose coupler curve interpolates nine general points in the plane, namely 1442 [55] : finding 
the maximal likelihood degree for many cases of matrices with rank constraints [45] and 
observing duality which was proven in [26]; a range of results in physics such as [40114111541155] : 
and numerically solving systems of nonlinear differential equations [32l 139] . 

We consider the equation ([T]) where r is the generic rank and the vfis are unknowns. 
Starting from one decomposition for T, we can move T(s) along a loop, for 0 < s < 1, such 
that T(0) = T( 1) = T. This consequently defines corresponding vectors Vj(s) which satisfy 


T(s) = 22v' 1 (s)®---®v' d (s). 
2=1 


The decompositions of T at the end values, s = 0 and s — 1, may be different. Since 
this process is computationally cheap, it can be repeated with random loops a considerable 
number of times and one can record all of the distinct decompositions found. Moreover, 
in the perfect case, where decompositions correspond to solutions to system of polynomial 
equations with the same number of variables, i.e., a square system, one can use a-theory 
via alphaCertif ied [42,33] to prove lower bounds on the total number of decompositions. 
Theory guarantees that all decompositions can be found using finitely many loops while ex¬ 
perience shows that all decompositions of T can be found after a certain number of attempts 
using random loops. When the number of decompositions is small, this process stabilizes 
quickly to the total number of decompositions. We describe using this process on some pre¬ 
viously known cases and predict several unknown cases. In particular, the values reported 
in Section [3] are provable lower bounds that we expect are sharp. 
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To put these results in perspective, we recall that finding equations that detect tensors 
of small rank is a difficult subject. Recent progress is described in [3], which gives a semi- 
algebraic description of tensors of format (n,n,n) of rank n and multilinear rank (n,n,n). 
In addition, several recent algorithms and techniques are available to find best rank-one 
approximations m or even to decompose a tensor of small rank [I2|]14j[58]|6?]. However, 
the problem generally becomes more difficult as the rank increases so that decomposing a 
tensor which has the generic rank is often the hardest case. 

The formats (3,4, 5) and (2, 2, 2, 3) were exceptional in our series of experiments since our 
technique showed that they have a unique decomposition (up to reordering). Indeed, an 
adaptation of the approach developed in [5<5] allowed to us to confirm our computations. 

Theorem 1.1. A general tensor of format (3,4,5) has a unique decomposition (QJj as a sum 
of 6 decomposable summands. 

Theorem 1.2. A general tensor of format (2, 2, 2, 3) has a unique decomposition (QJj as a 
sum of 4 decomposable summands. 

These theorems are proved in Section 0 The proofs provide algorithms for computing 
the unique decomposition, which we have implemented in Macaulay2 [31]. Based on the 
evidence described throughout, we formulate the following conjecture. 

Conjecture 1.3. The only perfect formats (ni,...,rid), i.e., R(ni,...,nd) in (J2]) is an 
integer, where a general tensor has a unique decomposition HD are: 

(1) (2, k, k) for some k — matrix pencils, known classically by Kronecker normal form, 

(2) (3,4, 5), and 

(3) (2, 2, 2, 3). 

We would like to contrast the tensor case to the symmetric tensor case, where the excep¬ 
tional cases were known since the 19 th century, as well as the following recent result. 

Theorem 1.4 (Galuppi-Mella [2D])- The only perfect formats ( n,d ), i.e., n" 1 • ( n+ 1~ 1 ) is 
an integer, where a general tensor in Sym d C n has a unique decomposition are: 

(1) (2, 2k + 1) for some k — odd degree binary forms, known to Sylvester, 

(2) (3,5) — Quintic Plane Curves (Hilbert, Richmond, Palatini), and 

(3) (4,3) — Cubic Surfaces (Sylvester Pentahedral Theorem). 

Theorem II .41 was still stated as a conjecture (following [53] and [52]) in the first preprint of 
the present article. See [53] and [60] for classical references. In [58] § 4.4], two of the authors 
showed that the Koszul flattening method predicts exactly the cases listed in Theorem 11.41 
and no others. 

2. Some known results on the number of tensor decompositions 

2.1. General tensors. The following summarizes some known results about tensors of for¬ 
mat (ni,... ,Ud). For any values of r smaller than the generic rank, which was defined in 
the introduction, the (Zariski) closure of the set of tensors of rank r is an irreducible alge¬ 
braic variety. This variety is identified with the cone over the r th secant variety to the Segre 
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variety P(C ni ) x • • • x P(C nti ) of decomposable tensors, e.g., see [T61I49] . In particular, it is 
meaningful to speak about a general tensor of rank r. 

Throughout this section, we consider cases where d > 3 and, without loss of generality, 
assume that 2 < n\ < ri 2 < ... < n d - First, we review the known results on the so-called 
unbalanced formats. 

Theorem 2.1. For formats (n \,..., rid), suppose that rid > Ilf=i n i ~ Ef=] l i n i ~ !)■ 

(1) The generic rank is min (rid, nf=i n *)- 

(2) A general tensor of rank r has a unique decomposition if r < nf=i — EtiK- 1 )- 

(3) A general tensor of rank r = Y[iZi n i ~ Ef=i( n * — 1) has exactly (^) different de¬ 
compositions where 

D (n, - 1)! ■ • • (nj_! - 1)!' 

This value of r coincides with the generic rank in the perfect case: when r = n d . 

(4) Ifn d > nti n i ~ Etafdh - 1), « general tensor of rank r > fltfl n t - Ef=i ( n i ~ 1); 
e.g., a general tensor of format (n\,... ,nd), has infinitely many decompositions. 

Proof. When n d > YliZi n i ~ Ef=i( n * — 1), Item [T] follows from pj8, Tlnn. 2.1(1-2)] (see 
also [13:, Prop. 8.2]). In the perfect case, i.e., n d = nf=i rii — Ef= i( n i ~ 1); Item [T] follows 
from P31 Prop. 8.3]. Items [2 and [3] follow from p3, Prop. 8.3, Cor. 8.4], When n d — 1 > 
nf=i Item H] follows from [H Lemma 4.1]. If n d — 1 = nf=i 

then nti* ri * = Ef=i( n * — 1)- Hence, 1 + Ef=i( n * — 1) cannot divide nf=i an d s ° the 
format cannot be perfect. □ 

The case (2 ,n,n), corresponding to pencils of square matrices, is the only case for which 
the binomial coefficient (° r ) in Theorem 12.1 ll3l) is equal to 1. The unique decomposition is a 
consequence of the canonical form for these pencils, found by Weierstrass and Kronecker [IB], 
For convenience, Table [T] lists some perfect cases coming from Theorem I2.1HHD . namely 
when n d = IlEi ~ Ei=i( n * — 1) with generic rank r = n d - 

After Theorem 12.11 the only open cases are when the balancedness condition is satisfied: 

d— 1 d— 1 

(3) n d < Y\rij - y^frij - 1). 

i— 1 i= 1 

A seminal identihability result for general tensors up to a certain rank is (SHI Cor. 3.7]. 
In [2U], based on weak defectiveness introduced in PJJ], there are techniques to check the 
number of decompositions of a general tensor of rank r, generalizing Kruskal’s result [48 ] . 

De Lathauwer’s condition (mentioned after equation (1.7) in [25]) can guarantee unique¬ 
ness up to rank 

d- 1 

y; rii — d + 1 
1=1 


( 4 ) 
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(n i,.. .,n d ) 

gen. rank 

# of decomp, of general tensor 

(2, n, n) 

n 

1 

(3,3,5) 

5 

6 

(3,4, 7) 

7 

120 

(3,5,9) 

9 

5005 

(3,6,11) 

11 

352716 

(4,4,10) 

10 

184756 

(2, 2, 2, 5) 

5 

6 

(2, 2, 3, 8) 

8 

495 


TABLE 1. Generic ranks and numbers of decompositions for perfect formats (Thm. |2T|) . 


In the 3x4x5 the bound ([4|) is 5. In the 2 x 2 x 2 x 3 the bound (J4J) is 3. Indeed De 
Lathauwer considers only “tall” tensors when the rank is no larger than the dimension of any 
mode, so his methods don’t apply to the perfect case. However if you ignore this assumption 
(rank R < n^) in j2S Theorem 2.5], the second inequality 

becomes 30 < 36, which is valid in the 3x4x5 case. So the previous bound just misses our 
results. 

For all formats such that Y\i=i n i — 15,000 which satisfy the inequality ([3]), a general 
tensor of rank r which is strictly smaller then the generic rank has a unique decomposition 
except for a list of well understood exceptions, e.g., see (20]. Thm. 1.1], These results support 
the belief that, other than some exceptions, a general tensor of rank r smaller then the generic 
rank has a unique decomposition. When r is the generic rank, since the techniques in [2U j 
cannot be applied, we apply numerical algebraic geometry to such cases in Section 0 

2.2. The symmetric case. The following summarizes results about symmetric tensors to 
contrast with the general case. Recall that symmetric tensors of format (n, d) are tensors 
T e Sym d C n , which can be identified with homogeneous polynomials of degree d in n 
variables. The (symmetric) rank of T is the minimum r such that there is an expression 


T = 


v l <g) • ■ ■ (8) v l 


i— 1 


with v l G C n . If T is identified with a polynomial, then each summand v l <g) • • • <g) v l is the 
d-power of a linear form. By a naive dimension count, a general tensor in SynT'C” has rank 
at least n _1 ■ (" + ^ -1 )- When this fraction is an integer, the symmetric format (d, n) is called 
perfect. As in the general case, perfectness is a necessary condition in order for a general 
tensor in Sym d C n to have only finitely many decompositions. 

The following is the basic result about decomposition of symmetric tensors that we state 
for perfect formats. 
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d 

gen. rank 

ff of decomp, of general tensor 

4 

6 

oo 

5 

7 

1 

7 

12 

5 

8 

15 

16 


TABLE 2. Generic ranks and numbers of decompositions for general tensors 
in Sym d C 3 [60]. 


Theorem 2.2 (Alexander-Hirschowitz [3]). Let d > 3 and assume r = n 1 • ( n+ ^ G Z. 
A general tensor in Sym d C n has finitely many decompositions of rank r except in the follow¬ 
ing cases: (n,d) = (3,4), (5,3), or (5,4). In these three exceptional cases, a general tensor 
has no decomposition of rank r, but infinitely many decompositions of rank r + 1. 

When n = 2, note that G Z exactly when d is odd. In these cases, Sylvester proved 
that there is a unique decomposition with summands [6Uj. 

Theorem II.4f [50]) lists two other cases when a general tensor in Sym^CT has a unique 
decomposition, namely Sym 3 C 4 where (n,d) = (4,3) and Sym 5 C 3 where (n,d) = (3,5) [60] . 

When n = 3, note that | • ( a! ^ 2 ) G Z is an integer exactly when d = 1 or d = 2 modulo 3. 
Table [2] records all the cases that can be found in [60] concerning Sym d C 3 . The clever syzygy 
technique used in [60] seems not to extend to higher values of d. 

Remark 2.3. Let d = 1 or d = 2 modulo 3. By [23] Thm. 4.2(vi)], the number of decom¬ 
positions of a general symmetric tensor in Sym d C 3 is bounded below by the degree of the 
tangential projection from r — 1 points, where r = ((/+2 | ; ^ + b j s the g ener j c rank. This latter 
degree is computed as the residual intersection of two plane curves of degree d having r — 1 
double points, which is d 2 — 4(r — 1) = d 2 — 4 ^~ 1| | ! l/+ ' 11 = LAilLLT. 

An analysis of the degeneration performed in [25] suggests that actually the number of 
decompositions of general symmetric tensor in Sym d C 3 should be divisible by ( d ~ 2 W~ 4 ) . This 
guess agrees, for d < 8, with the above table from [60] and the results for d < 11 in 4‘i.31 

A general symmetric tensor of rank r which is strictly smaller then the generic rank has 
a unique decomposition except for a list of well understood exceptions, see [6] [2X1153], 

3. HOMOTOPY TECHNIQUES FOR TENSOR DECOMPOSITION 

In this section, we first describe the monodromy-based approach we use to determine 
the number of decompositions for a general tensor. The software Bert ini PMJ is then 
used in the subsequent subsections to compute decompositions for various formats. In the 
perfect cases under consideration, the number of decompositions can be certihably lower 
bounded via alphaCertif ied [421145] . In particular, for the two cases (3,4, 5) and (2, 2, 2, 3) 
which are discovered here to have a unique decomposition, we provide theoretical proofs 
in Section [5] Our computational methods include probabilistic reductions (e.g. cutting by 
random hyperplanes and choosing random points) and numerical computations which are 
always subject to round-off errors in any finite computation (e.g. numerical path tracking). 
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Even though these methods are now completely standard, have been carefully and repeatedly 
tested, and yield completely reproducible results, they are technically only true with high 
probability, or up to the numerical precision of the computers we use. Recently results whose 
proofs partially rely on such methods have been denoted Theorem* (see [59]). In this article 
we say that these results hold with high confidence. 


3.1. Decomposition via monodromy loops. In numerical algebraic geometry, mon- 
odromy loops have been used to decompose solution sets into irreducible components [H]. 
Here, we describe the use of monodromy loops for computing additional decompositions 
of a general tensor. For demonstration purposes, suppose that a general tensor of format 
(ni,..., nfi) has rank r and finitely many decompositions. 

The approach starts with a general tensor T of format (ni,... with a known decom¬ 
position (JT]) with Vj G C n '‘ for i — 1,..., r. In practice, one randomly selects the Vj first and 
then computes the corresponding T defined by (jT]) . To remove the trivial degrees of freedom, 
we assume that (v))i = 1 for i = 1,..., r and j = 1,..., d — 1. That is, we have a solution of 


Ft{v l-..,v r d ) = 


T ~ E[= i v\ ® • • 
(u?)i - !> 


®v l d 


i = 1 ,- 1 


= 0 . 


The system F T consists of Y]?j=i n i + r (d — 1) polynomials in r • Ej=i n j variables. Since 
r = R(ni ,..., nfi) in (l2|) . the number of polynomials is equal to the number of variables 
meaning that Ft is a square system. 

Now, suppose that S C (C ni x • • • x C nd ) r consists of the known decompositions of T. For 
a loop r : [0,1] —>■ C ni "' nd with r(0) = r(l) = T, consider the homotopy 

H(vl,...,v r d ,s) = F r{s) (vl,...,v r d ) = 0. 


The loop r is selected so that the solution paths starting at the points in S when s = 0 are 
nonsingular for s G [0,1]. This is the generic behavior for paths r since the singular locus 
is a complex codimension one condition while we are tracking along a real one-dimensional 
arc r(s) for 0 < s < 1. The endpoints, namely at s = 0 and s = 1, of these solution paths 
form a decomposition of T. If a new decomposition is found, it is added to S. The process 
is repeated for many loops r. We leave many details about path tracking to USES]. 

Since F T and the homotopy FI is naturally invariant under the action of the symmetric 
group on r elements, we only need to track one path starting from one point from each 
orbit. Each loop is usually computationally inexpensive so we can repeat this computation 
many times. Experience has shown that randomly selected loops are typically successful 
at generating the requisite monodromy action needed to obtain all decompositions starting 
from a single one in a relatively small number of loops. 

In the subsequent subsections, when an exact value is reported, this means that at least 50 
additional randomly selected loops failed to yield any new decompositions. Thus, we expect 
that these values are sharp. When a lower bound is reported, this means that we have 
terminated the computation with the last loop generating many new decompositions. Thus, 
these lower bounds are probably quite far from being sharp, but do show nonuniqueness. 
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(ni,ra 2 ,n 3 ) 

gen. rank 

ff of decomp, of general tensor 

(3,4,5) 

6 

1 

(3,6,7) 

9 

38 

(4,4,6) 

8 

62 

(4,5,7) 

10 

> 222,556 


TABLE 3. Results of our numerical computation for all perfect format 3- 
tensors satisfying ((3]) with nLi — 150 and the number of decompositions 
for the generic tensor. 


(ni, ...,n d ) 

gen. rank 

ff of decomp, of general tensor 

(2, 2, 2, 3) 

4 

1 

(2, 2, 3,4) 

6 

4 

(2, 2, 4, 5) 

8 

68 

(2, 3, 3,4) 

8 

471 

(2, 3, 3, 5) 

9 

7225 

(3, 3, 3, 3) 

9 

20,596 

(2, 2, 2, 2,4) 

8 

447 

(2, 2, 2, 3, 3) 

9 

18,854 

(2, 2, 2, 2, 2, 3) 

12 

> 238,879 


TABLE 4. Results of our numerical computation for all perfect format tensors 
with d > 4 satisfying (j3]) with nf=i — 100.and the number of decompositions 
for the generic tensor. 


3.2. The number of decompositions for perfect format tensors. 

Computation 3.1. Tables[ 3] and[|] summarize the results of our numerical computations 
which determine (with high confidence) the generic ranks and numbers of decompositions 
for general tensors d > 3 satisfying (j3J) with EUi^ 100 - T(l ble El records our results for 
symmetric tensors. 

The generic rank is known to be equal to the expected one for formats (n, n, n ) [5Tj, which 
is not perfect for n > 3, and (2,..., 2) for at least k > 5 factors [IT], which is perfect if k + 1 
is a power of 2. A numerical check for k = 7 shows it is not identifiable. 

3.3. The number of decompositions for symmetric tensors. We highlight a few cases 
for computing the number of decompositions of symmetric tensors. 

In the cases of this table, Theorem 1.1 in [52] and the recent [30] imply that the number 
of decompositions of a general tensor is at least 2. 

For Sym d C 3 and d — 1 or d = 2 modulo 3, the expectation stated in Remark 12.31 is that 
the number of decompositions is divisible by ( d ~ 2 K d ~ 4 ) _ This is confirmed for d — 10 with 
320 = 20-16 and d= 11 with 2016 = 96 • 21. 
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Tensor space 

gen. rank 

# of decomp, of general tensor 

Sym 10 C 3 

22 

320 

Sym 11 C 3 

26 

2016 

Sym 5 C 4 

14 

101 

Sym 3 C 7 

12 

98 


Table 5. Results of our numerical computation for the number of decompo¬ 
sitions of generic tensors for some symmetric tensor formats. 


4. Pseudowitness sets and verification 

The approach discussed in Section [3] uses random monodromy loops to attempt to gen¬ 
erate new decompositions. Clearly, when showing that a format is not identifiable, one 
simply needs to generate some other decomposition. We can use the numerical approxima¬ 
tions to generate a proof that it is not identifiable in the perfect case using, for example, 
alphaCertified [121145] . However, to determine the precise number of decompositions, 
we simply run many monodromy loops and observe when the number of decompositions 
computed stabilize. In this section, we describe one approach for validating the number of 
decompositions and demonstrate this approach in Section 14.21 for counting the number of 
decompositions for a general tensor of format (3, 6, 6) of rank 8. 


4.1. Using pseudowitness sets. For demonstration purposes, consider counting the num¬ 
ber of decompositions of a general tensor of format (rii,..., rid) of rank r. Consider the 
following where we have removed the trivial degrees of freedom by selecting elements to be 1: 


G : = 



,v r d ) 


T = ELi v \ ® ■ • ■ ® v\, 

{vfii — 1 for i — 1 ,.. ., r and j — 1 ,..., d — 1 


The graph G is clearly an irreducible variety. Hence, the image n(G) is also irreducible 
where 7r(T, v \,..., v \,..., v [,..., v%) = T. If dim G = dim 7r(G), then we know that a 
general tensor T of format (ri\,... ,rid ) of rank r has finitely many decompositions, namely 

|7r~ 1 (r) n g\ 

t\ 

In particular, |7 t _ 1 (T) flG| is the degree of a general fiber of 7r with respect to G and the 
denominator r! accounts for the natural action of the symmetric group on r elements. 

Using numerical algebraic geometry, computations on G will be performed using a witness 
set , e.g., see [lflU66j . and on 7r(G) using a pseudowitness set [46U47j . A byproduct of computing 
a pseudowitness set for 7r(G) is the degree of a general fiber. 

Suppose that X C CA is an irreducible variety of dimension k and / is a system of 
polynomials in N variables such that X is an irreducible component of the solution set 
defined by / = 0. Then, a witness set for X is the triple {/, £,U} where C C C N is a 
general linear space of codimension k and V — X D C. Here, “general linear space” means 
that C intersects the smooth points of the reduction of X transversely so that \V\ = degX. 
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To focus only on the case of interest, we assume that 7r : C N —> C m is the projection map 
onto the first m coordinates such that 7 r is generically Gto-one on X , i.e., dim X = dim 7r(X). 
Then, a pseudowitness set for 7r(X) is the quadruple {f,7r,Xi,W} where W = X D M. 
and M. = x C N ~ m with C C m being a general linear space of codimension k. 
Here, “general linear space” means that Xi n intersects 7r(X) transversely, so that |7r(kF)| = 
|7r(X) fl M w \ = deg7r(X), and, for each T G 7r(X) D X4 w , |7t _ 1 (T) fl X\ — i. Therefore, n is 
an Gto-one map on W with \W\ = £ • deg7r(X). 

Returning to the problem at hand, since we aim to compute a pseudowitness set for 7 r(G), 
we can simplify this computation by only considering the fiber over a general curve section 
of 7 t(G). If k = dim 11 (G), let A C C ni ® ■ ■ ■ <S> C nd be a general linear space of codimen¬ 
sion k — 1. Then, G = 7r _1 (7r(G) flM) flG and 7 t(G) = 7t(G) C\A are both irreducible curves, 
e.g., see [B21 Thm. 3.42] and [SHI Thm. 13.2.1], respectively. Moreover, deg7r(G) = deg7r(G) 
and the degree of a general fiber for G and G with respect to 7r are equal. 

With this, we now aim to compute a pseudowitness set for the curve 7 t(G). Since G is a 
curve, we compute a pseudowitness set for 7 t(G) by intersecting G with a hyperplane of the 
form 7 r" 1 ('H) where 1~L C C ni <E> • • • <g) C nd is a general hyperplane. Since 7 r _1 ("H) is invariant 
under the natural action of the symmetric group on r elements, we will first consider the 
intersection of G with a general element of an irreducible family of hypersurfaces that are 
invariant under the same action and contains hyperplanes of the form 7 r” 1 ('H). This is 
sufficient by the results of [SB] and simplifies the computation since only one point in each 
orbit needs to be computed, i.e., only one point out of every r! based on the natural action 
of the symmetry group on r elements. 

Remark 4.1. At first sight, by reducing down to the curve case, it might seem that some 
technicalities could be avoided by choosing a general hyperplane among the hyperplanes 
invariant under the symmetric group and only needing to compute one point in each orbit. 
However, there could be some difficulties that arise with this. 

The first difficulty is that the set of hyperplanes invariant under a finite group does not 
have to be irreducible. For example, let GcC* denote the sixth roots of unity. Consider 
the action of G on C 2 defined by g ■ ( z,w ) ha (g 2 z,g 3 w). In this case, the set of invariant 
hyperplanes consists of two points, i.e., the 2 axis and the w axis. 

Even when there is a large family of invariant hyperplanes that make sense, the second 
difficulty is that none of them need be invariant enough to intersect the algebraic set of 
interest transversely in the degree number of points. For example, let X C C 2 be the 
solution set of z 2 — w 3 which is a curve of degree 3. Let G denote {1,-1} which acts on 
C 2 by g ■ ( z,w ) ha ( gz,w ). Note that X is invariant under G. The G-invariant hyperplanes 
consist of two components. One component is made up of the fibers of the projection map 
(z, tc) 1 —y w which meet X in two points. The other component consists of the w axis which 
is not transversal to X. 

To avoid these potential difficulties, we first define an irreducible family of hyperplanes that 
contains the invariant hyperplanes of interest. Then, the general theory in [5B| shows that 
this process computes the requisite points needed in the construction of a pseudowitness set. 
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4.2. Tensors of format (3, 6, 6) of rank 8. The tensors of format (3, 6, 6) have generic 
rank 9 in which a general tensor of this format has infinitely many decompositions. In [22], 
the open problem of computing the number of tensor decompositions of a general tensor of 
rank 8 of format (3, 6, 6) was formulated. To the best of our knowledge, this is probably 
the last open case when a generic tensor of some rank strictly smaller than the generic one 
is not identifiable. Theorem 3.5 of [22] proved that the number of decompositions is > 6. 
Moreover, [22] showed that the number of decompositions of format (3, 6, 6) of rank 8 is equal 
to the number of decompositions of a general tensor in Sym 3 C 3 ® € 2 ® C 2 , which is perfect 
with generic rank 8. We use the approach from £14.11 to show exactly 6 decompositions. 
Another approach based on a multihomogeneous trace test in [04] confirms this result. 

To use the approach from £ 14.11 consider the Veronese embedding 


v 3 : C 3 —y C 10 where ( x , y, z ) H > (x 3 , x 2 y , x 2 z, xy 2 , xyz, xz 2 , y 3 , y 2 z, yz 2 , z 3 ). 


We picked a random line C C C 40 and consider the irreducible curve 


C -.= 


^ (-T Oi, bx , ci,..., og, 6g, eg) 


T = X]i = i v 3 (cii) <g> <8> Ci € C 

K)i = (6di = 1 for i = 1,..., 8 


C C x (C 3 x C 2 x C 2 ) 8 . 


To compute a pseudowitness set for we need to compute C D (Ain x (C 3 x C 2 x C 2 ) 8 ) 

where Ai n C C 40 is a general hyperplane. Consider the irreducible family of hyperplanes 
defined by the vanishing of linear equations of the form 

40 3 8 2 8 2 8 

(5) J2 a i T i + Pi ^2( b k)j + S i J2^i = e ‘ 

j =1 j =1 k =1 j =1 fc=l j =1 k =1 

By construction, this family is invariant under the natural action of the symmetric group Sg 
and contains all hyperplanes of the form Ain X (C 3 x C 2 x C 2 ) 8 . After picking a random 
hyperplane "H of the form ([5]) and starting with one point on CCXH, we used monodromy loops 
via Bertini to compute additional points in Cn'H which stabilized to 1020-8! = 41,126,400. 
The trace test [62] confirms that this set of points is indeed equal to CfVH. After selecting a 
random hyperplane Ai n , we then computed W = C D (M.n x (C 3 x C 2 x C 2 ) 8 ) by deforming 
from C nV. which yielded |W| = 6 • 8! = 241,920. Since n(C) D Ai n = jC which has degree 1, 
|7 r (W)| = 1 thereby showing exactly 6 decompositions. 

We summarize the result of this computation. 

Computation 4.2. Numerical algebraic geometry together with randomly selected linear 
spaces show (with high confidence) that a general tensor in Sym 3 C 3 (8) C 2 <g) C 2 has exactly 6 
decompositions. 


5. Tensor decomposition via apolarity 

In [5E], a technique generalizing Sylvester’s algorithm was implemented by considering the 
kernel of the catalecticant map, which in turn is a graded summand of the apolar ideal. In 
principle, this apolarity technique can be used for any algebraic variety. 
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5.1. A uniform view of (Koszul) flattenings. Let V and V, be arbitrary finite dimen¬ 
sional vector spaces over C of dimensions n and n*, respectively. For — n < p < n, let V 
denote fundamental representations of GL(V) where we interpret /\ p = A P K* when p < 0. 
For a multi-index / G 7L d , let Vj denote the tensor product of fundamental representations 

v, ■.= ... ®r( d v d . 

Note that V v i := ... ®V d . We may assume, up to reordering, that ij > 0 for j — 1,..., h, 

ij < 0 for j — h + 1,..., d. We obtain linear maps K p \ flfV —» fl( +1 V that depend linearly 
on V by way of the Koszul complex. Specifically, for v G V and <p G flfV define 

K p (v)(tp) = ip A v for p > 0, 

Kp(v)((p) = <p(v) for p < 0. 

Now we consider the tensor product of many Koszul maps, which are linear maps on tensor 
products of fundamental representations that depend linearly on T e := Vi<E>... ®V d . 

Kj(T): V!^V I+1 *. 

For indecomposable elements i>i<g)... (8 w d G ... ®V d and Lpi® ... ®ip d G Vj define 

h d 

(6) ... ®<p d ) = ($£)(<Pj A vj)® (££) (pj(vj)). 

3=1 j=h +1 

The definition of Kj is extended by bi-linearity. We often drop the argument (v) or (T) 
of K p or Kj when the tensor to be flattened is not specified and the linear dependence is 
understood. From this definition it is clear that the image of Kj{v\® ... ®v d ) is isomorphic to 

h d 

0 (A !j C Vj/{vj)) ® & (A“ v_1 ('«/)) • 

j =1 j=h +1 

A consequence of dimension counting, bi-linearity of Kj, and sub-additivity of matrix rank 
is the following, which is essentially described in [50], Prop. 4.1]. 


Proposition 5.1. Suppose T G has tensor rank r. Let ij > 0 for j — 1,..., h, ij < 0 

for j — h + 1,..., d. Then the Koszul flattening Kj(T ): Vj —> V I+1 d has rank at most 


h 

r i f ■ J 

3 =1 


Uj - 1 


d 

n 

j=h +1 


Uj ~ 1 

-ij - 1 


In particular, the ( r d + 1) x (77 + 1) minors of Kj(T) vanish. This is meaningful provided 
that r d < rriin{dim Vj, dim Vj + \d}. 


Let m = UU (A X ) ' U d j=h+i (-i— i)- The proposition says that for all tensors T with 
rank(T) = r we have rank(AT/(T)) > r ■ m. Thus, Koszul flattenings potentially provide the 
most useful information whenever the following ratio is maximized: 

min {dim V}, dim Vj +1 d} / m. 

When rm < min {dim Vj, dim V I+ id) we say that the flattening can detect rank r since we 
expect that the flattening KflT) will have different ranks when T has rank either r — 1 or r 
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(this is true in all the examples we have tried, including the ones occurring in Theorems 11.11 
and ll.2p and there is no obstruction to this based solely on the size of the flattening matrix. 

5.2. Apolarity Lemma for Koszul flattenings. Recall from (J6]) that 

h d 

Ki(vi® ... ®v d ) = 0 A Vj)® (g) ((pj{vj)) = 0 

1=1 j=h+ 1 

for all pure tensors (f G Vj. 

It is useful to look at tensors in the kernel of Kj(T) as linear maps. With this aim, we 
need to distinguish the negative and nonnegative parts of / G 7L d . So let N U P — {1,..., d} 
be the set partition such that —In G Zk( 0 , Ip G Z> 0 and the notation Ip (resp. In) is the 
vector in 7L d gotten by keeping the elements of / in the positions P (resp. N) and setting 
the rest of the entries to zero. We also let lp denote the vector with ones in the positions 
denoted by the index P (and zero elsewhere), and similarly for In- With this, we may 
identify Vj = Vj p+ j n = Hom(R_/ JV , V dp ), and consider the Koszul flattening of T G R(i,...,i) as 

K t (T): Rom(V_ lN ,V Ip )^Kom(V_ lN+lN ,V Ip+lp ). 

Up to reordering the factors, Kj(T) is defined on decomposable elements \&)j£i N w j) by 


K-! (ifi® ... ®v d ) (ip) 



W; G Hom(VW v , V Jp ). 


In our setting, [50], Prop. 5.4.1] yields the following lemma (see (J9]) for a concrete case). 
Since this technique refers to a vector bundle, it could be called “nonabelian” apolarity, in 
contrast with classical apolarity which refers to a line bundle (see [501 Ex. 5.1.2] and [58] § 4]). 


Lemma 5.2 (Apolarity Lemma). Suppose T = v f<g)... ®v^. 

ker A/(T) D {ip G Hom(U_ /jv , V Ip ) \ ip (V- In+1n A <g) jeN v j ) A (® jeP vf) = 0 for s = 1,..., r}. 

Proof. Pick ip in the space on the right hand side of the inclusion. Choosing any Wj G V_j N+ \ N 
for every j G In, we have 


k i (■ ■ -® v d ) w (& w j ) = XA (® k- A ) A 


,s=l 


Kjei N 


•5=1 V/G/jV 


kie/p 


and each summand vanishes by the assumption. 


□ 


5.3. The 3x4x5 case. Let us denote the three factors as A = C 3 , B = C 4 , C = C 5 . The 
following are all possible relevant non-redundant Koszul flattenings (up to transpose), which 
all depend linearly on A®B®C. 

K( 0 , _!,_!): ( B®C)* A, K { _ i i0 ,_i): (A®C)* B, iP ( _i,_i, 0) : (A®B)* C 
K m -i) : C*®A -G B®fiA , A ( _ 1 , 0 , 2 ): A*®fiC —>■ B®f\C 
A(o,i,-i) : C*®B —>■ A®/\ B , A(_ 1 , 1 , 0 ) : A*®B — y C®f\ B, 
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AT ( —!,o,i) = A*®C -A- B®ftC , A ( o,-i,i): -A A® AA 

Consider 

K U (T ): ft 1 A®ft 2 B®ft 3 C -A A u1+1 ^®A U2+1s ®A“ 3+1 C'. 

where, for any vector space V", we interpret negative exterior powers by asserting that 
AV = A"V* if s < 0. 

For example K 0 _i_i{a®b®c) has image 

{ftA A a)®{B*{b))®{C*{c)) C ftA®ftB®ftC. 


The factor {B*{b))®{C* (c)) is just any scalar, obtained by contracting b with B*, and c 
with C*. We are left with {ftA A a) = (a), which is 1-dimensional. 

Another example is K 0 i_i{a®b®c) which has image 

(A'U A a)®{ftB A b)®{C*{c)) C ftA®ftB®ftC. 

The factor C*{c) is a scalar obtained by contracting c with C*. This leaves {ftA A a) = (a) 
tensored with {ftB A b) C ftB, but {ftB A b) = {B/b)®(b), which is 3-dimensional. 

In general, the image of K u {a®b®c ) has dimension 


/dim A — 1\ /dim B — l\ /dim C — l\ 

V f{u i) /V f{u 2 ) /V f{u 3 ) )' 


where f(x) — I ^ X ^ . On the other hand, the maximum rank that K u can 

M ; [ -1 - 1 if 1 < 0 ’ u 

have is the minimum of the dimensions of the source and the target, or 


( 8 ) 


f f dim f dim B\ f dim C\ f dim A \ f dim B \ f dim C \ j 

mm iv l«il A l«2| )\ M J’Ui + ilJ U2 + 11JU3 + 11J;' 


Since the matrix rank of a Koszul flattening of a tensor T is bounded above by a multiplicative 
factor (J7J) of the tensor rank of T, the maximum tensor rank that a Koszul flattening can 
detect by a drop in matrix rank is the ratio of (JS]) and (J7J). For convenience we record the 
dimensions and the multiplication factor (jTJ) for each flattening. 


map 

size 

mult-factor 

max tensor rank detected 

A( 0 ,-i,-i) 

3 x 20 

1 

3 

iA-i,o,-b 

4 x 15 

1 

4 


5 x 12 

1 

5 


15 x 12 

2 

6 

K ( 1 , 0 - 1 ) 

12 x 15 

2 

6 

K( o,i,-i) 

18 x 20 

3 

6 

-ft (- 1 , 1 , 0 ) 

12 x 30 

3 

4 

-ft (- 1 , 0 , 1 ) 

40 x 15 

4 

4 

-ft ( 0 ,- 1 , 1 ) 

30 x 20 

4 

5 

-ft (-1,0,2) 

40 x 30 

6 

5 

-ft (0,-1,2) 

30 x 40 

6 

5 

-ft (0,-1,2) 

30 x 40 

6 

5 
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We see that the only maps that might distinguish between tensor rank 5 and 6 are Ayi,_i,o)) 
Ar ( i )0l —i), and Ay 0 .i,-i)- Since f\A = A*, the first two maps are transposes of each other: 

o) = (^d,o,-i))*- 

Thus, we proceed by considering Ay 1)0 _i) and Ay 0 ,i,-i)- 

In our case, Apolarity Lemma 15.21 says that 

(9) ker Ad,o,-i(Z)i=i aAc) D {y G Hom(C, A)\ip(cf) A a* = 0 for i = 1,..., s}. 
and 

ker A" 0 ,i,- 1 (^i= 1 a Ac) {p G Hom(C, B)\ip(ci) A b t = 0 for i = 1,..., s}. 

Equality should hold for honest decompositions, see [50] Prop. 5.4.1], 

With this setup, we now present a proof of Theorem 11.11 

Proof of Theorem U. 11 For general / G A(&B ®C , Ad,o-i(/) is surjective and ker Ad,o,-i(/) 
has dimension dim Hom(C, A) — dimA 2 A ® B = 15 — 12 = 3. To complete the proof 
we interpret the linear map Ad,o,-i(/) as a map between sections of vector bundles. Let 
X = P(A) x P(A) x P(C7), endowed with the three projections tia, tb, n c on the three 
factors. We denote 0(a,/3,^f) = n\0(a) <E> n^O(/3) <E> ttqO^). Let Qa be the pullback of 
the quotient bundle on P(A). 

Let E = Qa <S> 0(0, 0,1) and L = 0(1, 1,1). Note that E is a rank two bundle on X. As 
in [50l[58j, the map Ad,o_i(/) can be identified with the contraction 

Ab, 0 ,-i(/): H°(E) — >H°(E* ® L)* 
which depends linearly on / G H°(L)*. 

The general element in H°(E) vanishes on a codimension two subvariety of X which 
has the homology class 02 (E) G H*(X,1f). The ring H*(X,Z *) has three canonical gen¬ 
erators tA,ts,tc and it can be identified with Z[£a, A, tc\/(t% t%, t b c ). Since the rank of 
Qa is 2 and 0(0,0, 1) is a line bundle, we have the well known formula 02 (E) = C 2 (Qa) + 
ci(Qa)ci(O(0, 0,1)) + c\(0( 0, 0,1)). We have ci(Q a ) = t A , c 2 (Qa) = t\, c x (O(0, 0,1)) = t c 
and we compute 02 (E) = t\ + £Ac + Hence three general sections of H°(E) have 
their common base locus given by 02 (E ) 3 = (t 2 A + Ate + t^ff = This coefficient 6 

coincides with the generic rank and it is the key of the computation. Terracini’s lemma 
(see [49] Cor. 5.3.1.2]), using a simple tangent space computation, verifies that the generic 
rank is indeed 6. We pick a tensor / = aA<A constructed with six random points 

(a*, hi, of) G P(A) x P(A) x P(C) for % — 1,..., 6. By using Macaulay2 (see the Macaulay2 
hie attached at the arXiv version of this paper) we may compute that the variety 

{(a, c) G P(A) x P(C) A(c) A a = 0 \/cp G ker Ad,o,-i(/)} 

consists of the union of six points (ai,cf) for i = 1,...,6, and this holds even scheme- 
theoretically. It follows that equality holds in (jUJ) with s — 6, indeed Ad,o,-i(/) is surjec¬ 
tive and has 3-dimensional kernel, while the right-hand side of Q has dimension at least 
dim Hom(C, A) -6-2 = 3. Moreover the common base locus of ker Ad, 0,-1 (/) is given 
by 6 linear spaces {a*} x P(£>) x {q} for i — 1,..., 6. 

We claim that the common base locus of ker Ad,o,-i(/) is given by 6 linear spaces as above 
for general tensor /. Indeed, the common zero locus of three sections of E can be seen 
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as the zero locus of a section of E ® 3 . The dimension of the zero locus of a section of E ® 3 
is at least 9 — 6 = 3 and it is upper semicontinous with respect to the section, hence the 
common base locus of ker Ai,o,-i(/) is a pure-dimensional 3-fold for general /. Since the top 
Chern class of E 03 is (A) 3 , which we computed to be 6 t\t^, we know that the common 
base locus of ker Ad,o,-i(/) comes from a O-dimensional scheme of degree 6 on P(A) x P(C), 
after a pullback with 7 tb, by [29j Prop. 14.1 (b)]. Since this O-dimensional scheme consists 
of 6 distinct points for the tensor / that we have picked, this property remains true for 
general /. 

In particular, the decomposition / = Y^=i a* ® ® q has a unique solution (up to scale) 

for Oj and q. After a*, q have been determined, the remaining vectors fe* can be recovered 
uniquely by solving a linear system. □ 

Remark 5.3. If we attempt to repeat the same proof using Aq,i-i in place of AT 1)0 _i most 
parts go through unchanged. The map Ao,i_i(T): C*®B —* A®/\B is 18 x 20, and general 
element T produces a 2-dimensional kernel. Then, we consider the intersection of two general 
sections of H°(E) = Hom(C', B), where now E = Q B ®0( 0,0,1). The top Chern class of E 
is (by a similar calculation as in the proof of Theorem 11,11) 

4 t z B t z c + 3 t 2 B tQ. 

This gives that the common base locus of ker Ao,i_i(T) is given by a degree 7 curve on the 
Segre product Seg(P(C) x P (B)). This curve necessarily contains the 6 points needed to 
decompose T, but information from another Koszul flattening is needed to find them. 

Remark 5.4. For the (3,4, 5) format, we can even decompose a general tensor T of any rank r 
between 1 and 6. The trick is to add to T the sum of 6 — r general decomposable tensors, 
find the unique decomposition with the algorithm described in the proof of Theorem 11.11 
and subtract the 6 —r tensors that have to appear in the decomposition. Unfortunately, this 
technique cannot work in other cases if we do not have a tensor decomposition to start with. 

5.4. The > 4 factor case. We have seen that general tensors of format (2, n, n) and (3,4, 5) 
are identifiable. We asked if there are other formats with this property when there are > 4 
factors. To our surprise, the numerical homotopy method predicted an additional case where 
identihability holds. Our construction of Koszul flattenings and multi-factor apolarity above 
allows us to prove this fact, which we prove next. 

5.5. The 2 x 2 x 2 x 3 case. For this part, let A = B = C = C 2 and D = C 3 . Because of the 
small dimensions we are considering, the number of interesting Koszul flattenings for tensors 
in A®B®C®D is limited to the following maps, which depend linearly on A®B®C®D. 

The 1-flattenings (and their transposes): 

A.—i,o,o,o - ^4* —t B®C®D, Ao,—i,o,o: B* —> A®C®D, 

Ao,o,—i,o : C* —t A®B®D, Ao,o,o,-i ; D* —> A®B®C, 
which detect a maximum of rank 2 in the first 3 cases and a maximum of rank 3 in the last. 

The 2-flattenings (and their transposes): 

Aq,o,— i ,—i : C*®D* —>■ A®B, Ao,_i,o,-i: B*®D* —>■ A®C, 
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K— 1 , 0 , 0,—1 • A* ®E)* —> B®C. 

The maps are all 4x6 and detect a maximum of tensor rank 4. 

Remark 5.5. For (2,2,2,2), it is known that only two of the three 2-flattenings are algebraically 
independent, and the dependency of the third on the other two is “responsible” for the 
defectivity of the 3 rd secant variety <73 (P 1 x P 1 x P 1 x P 1 ). This secant variety has one 
dimension less than expected. This type of Segre variety was indeed studied by C. Segre [61j. 

The higher Koszul flattenings: 

K- 1 , 0 , 0 ,!: A*®D ->• B®C®f\D, AT 0 1 , 0 , 1 : B*®C ->■ A®C®f\D, 

-KojO,— 1,1 • C*®D —> A®B®f\D 

These maps are all 12 x 6 , and detect a maximum of rank 3. 

We will proceed with the 2-flattenings in the following proof of Theorem 11.21 since they 
are the only flattenings that detect the difference between rank 3 and 4. 

Proof of Theorem 11. HI Terracini’s lemma (see [49] Cor. 5.3.1.2]), using a tangent space com¬ 
putation, verifies the well known fact that the generic rank is 4. Suppose T G A®B®C®D 
is general among tensors of rank 4 and write T = ^,=1 a i®frj< 8 )Cj(g)dj. 

First consider the case iF 0 i o,-i,-i: C*®D* —y A®B. If T is general of rank 4, then 
Jl 0 ,o.— 1 ,—i(T) has rank 4, and must have a 2-dimensional kernel. Now, we apply (l9lh The 
points {ci®di} must be contained in the common base locus of the elements in the kernel 
of A'o i o,-i,-i(F). Let X = P(A) x P(F>) x P(C') x P (D), endowed with the four projec¬ 
tions tta, Tig, 7 Tc, 7T D on the four factors. We denote 0(a, fd, 7 , S) = 1 r\0(a) <S> ti%0(/3) ® 
<g) 7 T * d O{8). Consider the line bundle E = (9(0, 0,1,1) and L = 0( 1,1,1,1) over 
X. The ring H*(X,7j) has four canonical generators tA,ts,tc,tD and it can be identified 
with Z[tA,tB,tc,tD]/(t‘A,'t%,tc^D)- We have C\(E) = tc + to- As in the proof of Theo¬ 
rem o the common base locus of the 2 -dimensional kernel of iFo,o,-i,-iCO may be seen as 
the zero locus of a section of E ® 2 = 0(0, 0,1, l)® 2 , which has top Chern class c 2 (A). Since 
c\(E ) = 2 tctc + tf, we get c\(E)tc = 1 -tctjj, c\(E)t£, = 2 -tct 2 D . It follows that two general 
sections of E have common base locus given by a cubic curve, denoted Cq,d of bi-degree 
(1, 2) on Seg(P(C) x P(ZD)), pulled back to a codimension 2 subvariety of X. The projection 
to P(D) is a conic, which we denote Qc- 

Similarly for the next 2-flattening, K 0 - ll 0 ,-i : B*®D* —> A®C, we repeat the same pro¬ 
cess, where all the dimensions and bundles are the same except for a change of roles of C 
and B. By the same method we obtain another conic Qb in P (D). 

Finally, if Qc and Qb are general, Bezout’s theorem implies that they intersect in 4 
points in P(-D), say {[hi], [^ 2 ], [^ 3 ], [dj]}. Like in the proof of Theorem 11.11 we can check, 
using a Macaulay2 script, that starting from T = J2t=i a A c *dj with (aj, bi, Ci, df) G X 
for i = 1,... ,4, we get that the conics found by the above procedure intersect exactly in 
{[di], [d 2 ], [d 3 ], [d 4 ]}. By semicontinuity, like in the proof of Theorem ll.il . we have that the 
intersection between the common base locus of the kernel of iF 0 > o,-i,-i(T') and the common 
base locus of the kernel of /l 0 ,-i,o,-i(d”), corresponding to the class (2tctD+t 2 D )(2t B t£)+t 2 D ) = 
At B tct 2 D , consists of the pullback under 7 ta of four distinct points on ¥(B) x P(C) x P(ZD), 
for general T, namely of the 4 linear spaces P(A) x {bi} <g) {q} x {d*}. In particular, the 
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decomposition T = JA =1 a* <E> bi ®Ci®di has a unique solution (up to scale) for bi, Ci and dp 
After bi, Ci , di have been determined, the remaining vectors a* can be recovered uniquely by 
solving the linear system T = J2i=i a‘i®b i ®c i ®d i . □ 

6. Conclusion 

By using a numerical algebraic geometric approach based on monodromy loops, we are 
able to determine the number of decompositions of a general tensor. Since this approach 
determined that general tensors of format (3,4,5) and (2, 2, 2, 3) have a unique decomposi¬ 
tion, we developed explicit proofs for these two special cases. With the classically known 
generically identifiable case of matrix pencils, i.e., format (2 ,n,n), we conjecture these are 
the only cases for which a general tensor has a unique decomposition. 

We are currently researching other applications of this monodromy-based approach, in¬ 
cluding determining identifiability in biological models [8]. 
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